AD-A237  143 


NASA 

Technical  Memorandum  104361 


AVSCOM 
Technical  Report  90- C- 024 


Maximum  Life  Spur  Gear  Design 


M.  Savage  and  MJ.  Mackulin 
The  University  of  Akron 
Akron,  Ohio 


H.H.  Coe  and  J  J.  Coy 
Lewis  Research  Center 
Cleveland,  Ohio 


Prepared  for  the 

27th  Joint  Propulsion  Conference 

cosponsored  by  the  AIAA,  SAE,  ASME,  and  ASEE 

Sacramento,  California,  June  24-27, 1991 

■■c  riuiiic  : 

t'Qitev  j 

(\i/\sn 


us  ARMY 
AVIATION 
SYSTEMS  COMMAND 
AVIATION  UT  ACTIVITY 


91-02342 


91  6  17  089 


MAXIMUM  LIFE  SPUR  GEAR  DESIGN 


M.  Savage  and  B.a.  Mackulin 
The  University  of  Akron 
Akron,  Ohio  44325 


H.H.  Coe  and  J.J.  Coy 

National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 


A',- 

f  ' 

D'"' 


Tar 


!'  V  .•< 


a  i  i.  J*  i  C  •  1  1  C  L'a^.,  _ _ 


/ 


B/.  - . . 

Ci  1.T, rl  L'J  l  1  O'.i/ 

jAvall  &ai5/er 


Dl3  i 


Abstract 

Optimization  procedures  allow  one  to  design  a 
spur  gear  reduction  for  maximum  life  and  other  end 
use  criteria.  A  modified  feasible  directions 
search  algorithm  permits  a  wide  variety  of  inequal¬ 
ity  constraints  and  exact  design  requirements  to  be 
met  with  low  sensitivity  to  initial  guess  values. 
The  optimization  algorithm  is  described  and  the 
models  for  gear  life  and  performance  are  presented. 
The  algorithm  is  compact  and  has  been  programmed 
for  execution  on  a  desk  top  computer.  Two  examples 
are  presented  to  illustrate  the  method  and  its 
applicat ion. 

Nomenclature 

A  velocity  factor 

P  material  surface  constant,  psi 

b  design  parameter  scaling  constant 

C  design  problem  constant 

C  dynamic  capacity,  lb 

d  design  parameter  scaling  coefficient 

E  elastic  modulus,  psi 

F  force,  lb 

/  face  width,  in. 

Vf  feasible  direction  gradient  vector 

Vh  total  violated  constraint  gradient  vector 

J  AGMA  bending  stress  geometry  factor 
6 

I  life  (10  cycles  or  hours) 

M  merit  function 

Vm  merit  function  gradient  vector 

N  number  of  teeth 

R  gear  radius,  in. 

R  reliability 

RMS  surface  finish 

allowable  surface  strength,  psi 
As  optimizing  step  size 

T  temperature,  *F 

V  constraint 

V  velocity,  in. /sec 

V  stress  volume,  in.^ 

Vv  constraint  gradient  vector 

H  weight,  lb 

X  unsealed  design  vector 


X  AGMA  factor 

Y  scaled  design  vector 

z  depth  to  critical  stress,  in.  | 

a  angle  of  approach,  deg 

8  characteristic  life  (10®  cycles  or  hours) 

H  friction  coefficient 

1/  Poisson's  ratio 

p  surface  radius  of  curvature,  in. 

a  contact  stress,  psi 

shear  stress,  psi 
^  pressure  angle,  deg 

W  angular  velocity,  rad/sec 

Exponents 

b  Weibull  slope 

c  proportionality  factor 

h  proportionality  factor 

p  load-life  exponent 

Subscripts 

10  90  percent  probability  of  survival 

(reliability) 

ag  gear  addendum 

B  base 

b  bending 

c  center  distance 

d  total  dynamic 

f  flash 

g  gear 

H  Hertzian 

i  independent  parameter  index 

j  search  step  index 

k  constraint  index 

Li  parameter  lower  estimate 

pi  pinion  radius  to  gear  addendum  on  line  of 
action 

a  sliding 

t  tooth 

Ui  parameter  upper  estimate 

w  weight 


Introduct ion 


Tha  optimal  design  of  a  apur  gear  mesh  is  a 
problem  of  considerable  interest  in  mechanical 
design. Whan  fracture  of  the  gear  teeth  due  to 
bending  is  the  primary  mode  of  failure,  the  minimum 
number  of  teeth  which  avoids  interference  offers 
the  strongest  gear  set  for  a  given  size.’  However, 
as  speeds  increase,  so  do  the  prospects  for  pitting 
and  scoring  modes  of  failure.  Fitting  at  and  below 
the  pitch  point  on  the  pinion  tooth  is  a  failure 
mode  which  limits  the  life  of  the  gear  teeth.® 
Procedures  have  been  presented  to  design  a  gear  set 
for  minimum  size  and  minimum  weight,  considering^ 
pitting  fatigue  as  well  as  bending  and  scoring.  ' 
These  procedures  determine  an  optimal  gear  set  from 
a  weight  or  size  standpoint  alone. 

One  promise  of  optimization  is  adaptability  to 
the  application  at  hand.  For  example,  the  shaft 
center  distance  may  be  fixed  by  other  considera¬ 
tions  in  tha  design  of  a  machine.  Given  this  situ¬ 
ation,  one  would  like  to  apply  optimization  theory 
to  determine  a  design  which  maximizes  the  gear  set 
life  for  a  given  center  distance. 


as  specified  by  upper  and  lower  bounds  on  the  inde¬ 
pendent  design  parameters,  X,  such  that: 


-1.0  <  <  1.0 


(1) 


as 

Thus: 

Yj-djXj+bj  (3) 

where 

d^  -  - - -  (4) 

Xul  -  Xli 


and 


Xui  ♦  ^Li 


(5) 


The  computer  is  useful  here,  and  the  accessi¬ 
bility  of  the  personal  computer  makes  this  use 
attractive.  A  modified  feasible  directions  search 
algorithm  is  applied  in  a  continuous  design  space 
in  a  form  which  is  memory  and  time  efficient. 

Using  a  fixed  search  step,  which  is  halved  only 
when  the  search  passes  a  solution,  simplifies  the 
computations.  The  method  is  an  improvement  of  the 
basic  gradient  search  method.^'®  It  checks  the 
gradients  in  the  inequality  constraints  as  well  as 
the  merit  function  to  calculate  a  feasible  search 
direction  which  improves  the  merit  function  while 
staying  within  the  acceptable  design  region.  After 
finding  the  continuous  optimum  design,  the  program 
allows  the  designer  to  check  alternate  designs 
which  may  be  more  practical.  These  designs  can 
have  parameter  values  which  obey  the  discrete 
parameter  limitations  and  are  close  to  the  contin¬ 
uous  optimal  design's  parameter  values. 

This  method  is  applied  to  the  problem  of 
designing  for  maximum  life  at  a  given  center  dis¬ 
tance  as  an  example.  Also,  it  is  applied  to  a  sim¬ 
ilar  problem  of  minimum  center  distance  with  a 
given  life  requirement  for  comparison  purposes. 


The  actual  design  variable,  X^,  can  always  be 
retrieved  from  the  scaled  variable,  Y^,  by: 


Gradient 

Central  to  the  method  is  the  gradient  calcula 
tion.  This  is  performed  with  small  perturbations 
in  the  design  variables  from  the  nominal  position. 
The  gradient  in  the  merit  function,  Vm,  is  calcu¬ 
lated  as: 

9h 

9m 


9y 


Optimization  Method 


As  with  most  optimization  techniques,  the  pro¬ 
cedure  begins  with  several  vectors.  These  vectors 
are  the  independent  design  variables,  X;  the 
inequality  constraints,  V;  the  parameters  of  the 
merit  function,  P;  and  the  constants  which  define 
the  specific  problem,  C.  An  optimization  solution 
is  tha  design  variable  values,  X,  which  minimize  or 
maximize  the  merit  function  value  while  maintaining 
all  constraint  values,  V,  inside  their  specified 
limits.  A  procedure  starts  with  a  guess  for  the 
design  variable,  X,  and  iterates  with  some  logical 
procedure  to  find  the  optimal  design  variabie. 


where, 

9h 

M{Yj,...,Yj  *&Y . Y„)  -M(Yj . Y^ . Y„) 

5x  (8) 


Scaling 

To  maintain  balance  among  the  independent 
design  parameters,  the  design  space  is  scaled  into 
a  continuous,  dimensionless  design  space.®  The 
scaled  design  parameters,  Y,  vary  from  -l.*^  to  1.0 


In  Eq.  (8),  M  is  the  merit  function,  Y^  is  a 
scaled  design  variable,  and  iY  is  tha  small 
change  which  is  made  in  each  Y  .  In  tha  program. 
Ay  is  sat  at  0.001  which  is  O.OS  percent  of  the 
full  range  of  a  scaled  design  parameter. 
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Th«  magnituda  of  tha  gradient  vector  is  given 

by: 


|Vm| 


(9) 


For  minimization,  the  direction  of  change  in  Y 
which  reduces  the  merit  function,  M,  at  the  great¬ 
est  rate  is  determined  by  the  unit  vector,  Vm: 

Vm  -  -  (10) 

JW\ 

For  maximization,  the  sign  in  Eq.  (10)  reverses. 

In  the  simple  gradient  method,  Eq.  (10) 
defines  the  direction  for  the  step  change  in  the 
scaled  design  vector. 

-  Yj  +  As  Vm  (11) 

where  As  is  the  scalar  magnitude  of  the  step.  If 
no  constraints  are  violated,  this  will  be  the  next 
value  for  Y  in  the  search. 

step  size 

Step  size.  As,  is  a  significant  element  of  any 
optimization  procedure. ‘  For  stability  and  direct¬ 
ness,  a  fixed  step  size  is  used.  initially,  the 
step  size  is  0.1,  which  is  S  percent  of  the  range 
of  a  single  design  parameter.  But  the  procedure 
halves  the  step  whenever  a  local  minimum  is  reached 
or  the  search  is  trapped  in  a  constraint  corner. 

To  complete  the  simple  procedure,  the  search 
declares  a  solution  when  the  percent  change  in  the 
merit  function,  M,  is  less  than  a  pre-set  limit  of 
0.0001. 


<  0.0001 


(12) 


in  the  constraint  value,  V^,  for  motion  in  the 
gradient  direction.  The  vector  sum  of  the  gradi¬ 
ents  in  the  violated  constraints,  Vh,  is  the  second 
gradient  of  the  algorithm: 


k 

Ev-k 


(14) 


The  gradient  in  the  violated  constraints,  Vh, 
points  towards  the  acceptable  design  space  from  the 
unacceptable  design  space.  By  itself,  it  enables 
the  algorithm  to  turn  an  unacceptabie  initial  guess 
into  an  acceptable  initial  guess  by  a  succession  of 
steps: 

Yj.j  .  Yj  +  As  Vh  (15) 


Feasible  Direction 


Once  inside  the  acceptable  design  region,  the 
algorit)un  proceeds  along  the  steepest  descent 
direction  until  the  calculated  step  places  the  next 
trial  outside  the  acceptable  design  space.  To 
avoid  this  condition,  the  algorit)un  selects  a  fea¬ 
sible  direction  for  the  next  step.  Figure  1  shows 
a  sloped  constraint  intersecting  contour  lines  of 
improving  merit  function  values.  The  figure  shows 
gradients  in  the  merit  function,  Vm,  and  the 
impending  constraint,  Vh.  The  two  gradient  vectors 
ate  placed  at  the  last  viable  design  step,  although 
the  constraint  gradient  is  calculated  at  the  trial 
location  -  one  step  in  the  merit  gradient  direction 
from  the  indicated  point.  The  feasible  direction 
selected,  Vf,  is  the  unit  vector  sum  of  these  two 
gradients  in  the  merit  function,  Vm,  and  the  vio¬ 
lated  constraints,  Vh: 


Vf 


Vm  ♦  Vh 
|Vro  ♦  V(:j 


(16) 


Initial  Value 


And  the  next  step  becomes: 


The  optimization  procedure  described  above  is 
scaled,  fixed  step,  and  steepest  decent.  When  the 
initial  guess  is  in  the  acceptable  design  space, 
and  the  optimum  is  a  relative  minimum,  this  method 
wor)cs  quite  well.  However,  placing  the  initial 
guess  in  the  acceptable  r’.esign  region  is  often  a 
problem,  and  in  many  cases,  the  best  design  is 
determined  by  a  "trade-off"  among  conflicting 
design  constraints  at  the  edge  of  the  design  space. 


These  problems  are  addressed  with  a  second 
gradient.  Just  as  one  can  calculate  the  gradient 
in  the  merit  function,  one  can  calculate  the  gradi¬ 
ent  in  a  constraint  variable: 


(13) 


where  Vv  is  a  unit  vector  in  the  direction  of 
decreasing  value  in  the  constraint,  V^.  For  upper 
bound  constraints,  moving  through  the  design  space 
in  the  direction  of  Vv^^  will  reduce  the  con¬ 
straint  value  V^.  For  lower  bound  constraints,  a 
sign  reversal  in  Eq.  (13)  will  produce  an  increase 


Yj.j  -  Yj  ♦  As  Vf  (17) 


Alqorittun  Use 

By  using  subroutines  to  calculate  the  merit 
function  and  constraint  values  for  each  design 
trial,  the  procedure  separates  the  logic  of  the 
algorithm  from  the  analysis  necessary  to  define  the 
problem.  This  allows  the  design  problem  to  be 
changed  easily  without  concern  for  the  optimization 
procedure.  The  directness  of  the  procedure  adds 
additional  steps,  but  enables  the  program  to  run  on 
a  personal  computer. 

An  additional  benefit  of  separating  the  anal¬ 
ysis  routines  from  the  optimization  logic  is  the 
ability  to  modify  the  design  at  execution  and  ver¬ 
ify  the  characteristics  of  similar,  more  practical 
designs  with  the  same  program.  The  optimization 
procedure  wor)cs  with  a  continuous  design  space, 
which  includes  gears  with  fractions  of  teeth  and 
nonstandard  sizes.  By  allowing  the  user  to  see  the 
ideal  continuous  variable  solution  and  to  modify 
this  to  designs  with  whole  numbers  of  teeth  and 
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standard  sizes,  the  procedure  enables  a  designer  to 
determine  a  practical  optimum  design  easily. 

Input  Data 

The  algorithm  requires  the  user  to  provide  an 
ASCII  data  file  with  four  groups  of  data:  the  con¬ 
stants,  design  variables,  constraints,  and  merit 
function  weighting  coefficients.  All  constants, 
variables,  limits  and  coefficients  have  names  and 
units  in  the  data  file  to  make  the  program  output 
more  understandable.  In  addition,  the  data  file 
includes  a  title  for  the  optimization  task. 

The  constants  define  the  specific  problem  out¬ 
side  of  the  analysis  subroutines.  With  their 
labeling  names  and  units,  these  constants  define 
the  specific  problem  being  optimized  in  the  output 
data  file.  The  constants  also  enable  similar 
problems  to  be  solved  by  merely  changing  the  con¬ 
stant  values  without  recompiling  and  linking  the 
analysis  subroutines  to  the  optimizing  routines. 

The  design  variables  are  the  independent 
parameters  for  which  the  algorithm  seeks  values. 
Three  values  are  included  for  each  variable:  a 
lower  estimate,  an  upper  estimate,  and  an  initial 
guess.  The  lower  and  upper  estimates  determine  the 
scaling  range  as  the  span  between  them  is  set  to 
two  dimensionless  units  by  the  optimizing  proce¬ 
dure.  However,  they  do  not  set  hard  limits  for  the 
search.  If  hard  limits  are  needed  to  avoid 
indeterminate  calculations,  these  may  be  set  in  the 
constraint  limits. 


file  at  each  step  to  document  the  search  path  and 
assure  the  user  that  it  is  running  properly. 

Once  an  optimum  is  found,  the  program  writes 
the  design  variable  values  that  produced  the  opti¬ 
mum,  the  merit  function  value  and  its  component 
parameter  array,  and  the  full  constraint  variable 
array  to  the  screen  and  the  output  file.  It  pauses 
for  user  inspection,  and  on  resumption  gives  the 
user  a  chance  to  enter  a  different  set  of  design 
variable  values.  This  allows  the  user  to  include 
practical,  near  optimal  design  trials  in  the  com¬ 
puter  analysis  record.  Since  the  best  design 
should  be  of  this  type,  this  provides  full  documen¬ 
tation  to  the  additional  trials.  When  the  user 
declines  to  enter  mors  trials,  the  program  closes 
the  output  file  and  ends  execution. 


Analysis  Subroutines 

Two  analysis  subroutines  apply  the  optimi¬ 
zation  procedure  to  the  spur  gear  design  problem. 
The  routines  are  subroutine  BOUNDS  and  subroutine 
VALUES.  Subroutine  BOUNDS  analyzes  a  gear  design 
for  its  constraint  variable  values.  Subroutine 
VALUES  analyzes  a  gear  design  for  its  merit  func¬ 
tion  parameter  values.  These  subroutines  and  an 
input  data  file  treat  the  two  problems  of  designing 
for  maximum  life  at  a  given  center  distance  and  of 
designing  for  minimum  center  distance  at  a  given 
life. 

Bounds 


The  constraint  limits  include  a  bound  value 
and  a  direction:  lower  or  upper.  Both  a  lower  and 
an  upper  bound  may  be  set  on  the  same  property  by 
including  two  separate  limits  in  the  constraint 
list.  Since  the  program  reports  the  constraint 
values  for  the  optimal  and  check  designs,  one  can 
add  inactive  constraints  which  are  always  satisfid 
to  obtain  a  printout  of  additional  properties  in 
the  constraint  list.  This  also  gives  one  the 
opportunity  to  switch  the  controlling  constraints 
to  obtain  different  designs  for  different  physical 
conditions  without  changing  the  program. 

The  last  group  of  data  are  the  parameters 
which  may  be  included  in  the  merit  function  and 
their  weighting  coefficients.  These  weighting 
coefficients  may  convert  different  properties  to  a 
common  measure  such  as  cost.  They  may  be  order  of 
magnitude  corrections  to  keep  one  parameter  from 
dominating  the  design.  Or  they  may  all  be  zero 
except  one,  to  select  that  parameter  as  the 
quantity  to  be  optimized,  so  the  same  program  can 
be  used  to  obtain  different  designs  for  different 
requirements  by  switching  the  merit  function 
parameter  with  the  non-zero  coefficient.  This 
switch  changes  the  objective  of  optimization,  while 
the  switch  of  active  constraints  changes  the 
environment  of  optimization. 

Runtime  Options 

In  running  the  program,  the  user  is  requested 
to  enter  the  data  file  name  prefix.  The  program 
will  look  for  this  name  on  an  input  data  file  with 
the  suffix  ".IN"  and  will  write  the  output  to  a 
file  with  this  prefix  and  the  suffix  ".OUT"  as  well 
as  to  the  screen  as  it  runs.  Since  the  program 
calls  the  analysis  subroutines  many  times,  it 
writes  information  to  the  screen  and  the  output 


Included  in  the  16  constraint  variables  evalu¬ 
ated  by  subroutine  BOUNDS  are:  some  distances,  the 
pinion  weight  and  loading,  and  some  stress,  life 
and  scoring  limits.  The  involute  interference  var¬ 
iable  is  the  distance  along  the  line  of  action  from 
the  base  circle  of  the  pinion  to  the  addendum  cir¬ 
cle  of  the  gear.  This  distance  must  be  positive 
for  the  gear  tooth  tip  to  contact  the  pinion  tooth 
on  its  involute  surface  and  avoid  interference. 

The  pinion  weight  is  the  product  of  the 
density  of  steel  times  the  volume  of  a  disk  with  a 
diameter  equal  to  the  pinion  pitch  diameter  and  a 
thickness  equal  to  the  gear  and  pinion  face  width. 
The  transmitted  load  is  the  pinion  torque  divided 
by  the  base  radius  of  the  pinion,  which  is  the 
nominal  force  acting  between  the  gears  along  the 
line  of  action.  The  pitch  line  velocity  is  the 
rotational  speed  of  the  pinion  times  its  pitch 
radius. 

In  this  analysis,  the  dynamic  load  estimate  is 
the  AGMA  velocity  factor  model.  In  terms  of  a  gear 
quality  number,  Q.^,  the  AGMA  estimate  of  the  sum  of 
the  transmitted  load  and  the  dynamic  load  is: 


♦  /v 


(18) 


where 


A  -  50  ♦  56  1  - 


[12  -  QJ' 


4 


(19) 
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In  Bq.  (18),  ia  tha  total  dynamic  load,  is 
tho  nominal  transmittad  load  and  V  is  tha  pitch 
lina  valocity  of  tha  gaars.  In  Eq.  (19),  tha  qaar 
quality  numbar,  Q^,  may  hava  a  valua  batwaan  6  and 
11  with  11  corraspondinq  to  tha  hiqhar  quality 
qaar.  All  qaar  strassas  and  livas  are  calculated 
using  this  total  dynamic  load,  with  a  quality 
numtjar,  Q  >9. 


tha  load  sharing  factor,  is  a  tharmal  -  slastic 
factor  and  is  a  qaomatry  factor. 

Tha  pinion  lifa  and  mash  lifa  calculations  ara 
described  in  the  gear  life  model  section  which 
follows. 

Values 


Gear  tooth  banding  fatigue  and  gear  tip  scor¬ 
ing  ara  modes  of  failure  considered  in  the  program. 
The  banding  fatigue  model  uses  the  AGHA  J  factor  to 
estimate  the  banding  strass  with  tha  load  at  tha 
highest  point  of  single  tooth  loading  on  tha  pin¬ 
ion.  The  formula  for  the  bending  stress  is: 


(20) 


where  J  ia  the  AGHA  J  factor^*  and  f  is  the 
effective  gear  tooth  face  width. 

The  maximum  contact  stress  and  gear  tip 
Hertzian  pressure  are  calculated  as: 


F  ^  COS  ^ 

1  -  K  \/v\ 

“  + _ 1 

IP  / 

The  second  analysis  subroutine,  VALUES,  calcu¬ 
lates  the  three  components  of  the  merit  function: 
the  mesh  life  in  thousand  hours,  the  canter  dis¬ 
tance  in  inches,  and  the  pinion  weight  in  pounds. 
The  total  merit  function  is: 

M  -  C,»  ♦  C^(Rp  ♦  Rg)  ♦  C„W  (24) 

where  t  ia  the  mesh  life,  the  radii  sum  ia  the 
center  distance  and  W  is  the  pinion  weight.  The 
three  weighting  function  coefficients:  and 

are  constants  defined  in  the  input  data  file. 

By  assigning  different  values  to  these  coeffi¬ 
cients,  the  user  can  change  the  problem  being 
solved.  With  not  equal  to  zero  and  the  other 

two  coefficients  equal  to  zero,  the  optimization 
will  see)c  out  a  solution  which  maximizes  the  gear 
mesh  life.  With  and  equal  to  zero  and 

C  negative,  the  optimization  will  see)c  out  a 
solution  which  minimizes  the  canter  distance. 


Gear  Life  Model 


where  tha  p's  are  tha  radii  of  curvature  of  the 
two  gear  tooth  surfaces  at  the  point  of  contact, 
the  t/'s  are  the  Poisson  ratios  for  the  two  gear's 
materials  and  the  E's  are  the  elastic  moduli  for 
the  two  gear's  materials.  The  maximum  contact 
stress  occurs  at  the  lowest  point  of  full  load  con¬ 
tact  on  the  pinion  tooth.  The  gear  tip  Hertzian 
pressure  uses  one  half  of  tha  total  dynamic  load 
since  the  load  is  shared  between  two  tooth  pairs  at 
this  point. 

The  gear  tip  scoring  model  includes  the  pres¬ 
sure  times  velocity  factor  and  the  critical  oil 
scoring  temperature  model  from  lubrication  theory. 
The  normal  preasure  times  sliding  velocity  is  pro¬ 
portional  to  the  frictional  power  loss  of  the  gear 
set.  This  factor  ia  the  highest  for  contact  at  the 
gear  tip,  so  the  normal  pressure  is  the  gear  tip 
Hertzian  pressure,  the  sliding  velocity  at  the 
gear  tip  is  given  by: 

-  "p^pi  ♦  “pi> 


The  gear  life  model,  based  on  surface  pitting, 
comes  from  rolling  element  bearings.^  Surface  pit¬ 
ting  of  gear  teeth  follow  a  similar  pattern  to 
bearing  race  pitting,  with  the  possible  difference 
of  surface  initiation.  Lundberg  and  Palmgren  pro¬ 
posed  the  model  in  the  late  1930 's.^°  They  assumed 
that  the  log  of  the  reciprocal  of  the  reliability, 
R,  of  a  bearing  is  proportional  to  its  life,  ( ,  and 
some  stress  parameters.  These  parameters  are:  the 
stress  level,  T  ;  the  depth  to  the  maximum  shear 
stress,  z^f  and  the  stress  volume,  V.  The  rela¬ 
tionship  is: 


-h 

o 


Vf 


(25) 


In  relation  (25)  b  is  the  Heibull  slope  and 
c  and  h  are  exponents  of  proportionality  to  be 
found  experimentally.  This  is  the  equation  for  the 
two-parameter  Heibull  distribution  with  the  addi¬ 
tion  of  stress  and  size  factors.  In  general  terms, 
the  two-parameter  Heibull  distribution  is: 


where  the  W's  are  the  angular  velocities  of  the 
gears,  the  R's  are  the  radii  to  the  contact  point 
on  the  two  gears,  ^  is  the  nominal  pressure  angle 
of  the  gear  set  and  the  d's  are  the  angles  of 
approach  on  tha  two  gears. 

Tha  lubricating  oil  flash  temperature  la 
another  factor  used  to  monitor  gear  tooth  scoring. 
One  estimate  of  this  temperature^*  ia  given  by: 


(26) 


where  b  is  the  Heibull  slope  or  shape  factor  and 
0  is  the  characteristic  life  of  the  distribution. 
To  replace  the  characteristic  lifa  with  a  90  per¬ 
cent  probability  of  survival  life,  t solve 
Eq.  (26)  for  the  characteristic  life,  8. 


/I.x^a  l/v 

(23) 

0  -  1 

1 

i"rj 

(Rp  ♦  R,) 

Hi 

(27) 


where  T^  la  the  base  temperature  of  the  oil, 
is  the  average  surface  friction  coefficient,  Xp  is 
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At  R  •  0.9,  th«  Ufa  is  Substituting 

this  into  Bq.  (26)  givsa  tha  two  paramstsr  Wsibull 
distribution  in  tarms  of  tha  I  lifa  as: 


Lnfi]  .  LiI-L]  •  f-lLf 

"[r]  ■  lo.sj  ■  [»  J 


(28) 


Tha  lifa  to  raiiability  rslationship  of 
Eq.  (28)  is  for  a  spacific  load  which  determlnas 

the  I  lifa.  This  load,  F,  is  related  to  tha 

10 

component  dynamic  capacity,  C,  as: 


*  10 


(29) 


The  dynamic  capacity  of  tha  whole  gear  is 
lower  than  that  of  a  single  tooth.  In  a  single 
pass,  fixed  axis  gear  sat,  each  rotation  of  tha 
gear  subjects  every  tooth  on  tha  gear  to  a  single 
load  cycle.  Tha  gear  falls  whan  any  single  tooth 
on  tha  gear  fails,  and  tha  fatigue  damaga  in  each 
tooth  accumulates  independently  of  the  damage  in 
the  other  teeth.  In  successive  coin  tosses,  the 
probability  of  a  specific  combined  event  is  the 
product  of  the  probabilities  of  each  coin  toss.  So 
too,  the  raiiability  of  the  gear,  R  ,  is  the 
product  of  the  reliabilities  of  aac^  tooth  in  the 
gear. 


(33) 


Here,  the  dynamic  capacity  of  the  component,  C,  is 
the  load  which  has  a  90  percent  reliability  life  of 
one  million  cycles,  the  load  on  the  component  is 
F,  and  the  power,  p,  is  the  load-life  exponent. 
Since  the  life  at  the  dynamic  capacity  is  one 
million  load  cycles  or  unity,  it  does  not  appear  as 
a  variable  in  the  aquation. 

Because  gear  tooth  life  behaves  similar  to 
bearing  life,  engineers  at  the  NASA  Lewis  Research 
Center  formulated  a  model  for  gear  tooth  life  simi¬ 
lar  to  the  bearing  life  model. Starting  with 
Eq.  (25)  Coy,  Towqsend  and  zaretsky  developed  a 
model  for  the  reliability  and  life  of  a  spur  gear. 
The  model  uses  both  the  two-parameter  Waibull  dis¬ 
tribution  of  Eq.  (28)  and  the  Palmgrsn  load-life 
relation  of  Eq.  (29).  With  statistically  repli¬ 
cated  data,  they  showed  that  these  models  predict 
gear  tooth  pitting. 

From  Eq.  (25),  they  determined  a  relationship 
for  the  dynamic  capacity,  C^,  of  a  spur  gear  tooth. 
Rounding  the  exponents  to  whole  numbers  gives  the 
dynamic  capacity  as  a  function  of  Buckingham’s 
load-stress  factor, 

In  Eq.  (30),  B  is  a  material  strength  which  has 
the  dimensional  units  of  stress,  the  effective  face 
width  of  the  tooth  is  /,  and  the  curvature  sum  at 
tha  failure  point  for  the  contacting  teeth  is 
El//).  The  curvature  sum  is: 

Ll/p  »  _L  +  —  (31) 

P,  P, 


Here,  p  is  tha  radius  of  curvature  of  the 
gear  tooth  surface  at  the  failure  point,  and  p 
is  tha  radius  of  curvature  of  the  pinion  tooth 
surface  at  the  failure  point.  With  the  dynamic 
capacity  expressed  in  this  form,  tha  material 
strength  factor  serves  the  role  of  the  surface 
fatigue  strength,  of  the  AGHA  design  coda.‘* 

A  relation  for  tha  material  strength  factor  in 
tarms  of  the  surface  fatigue  strength  is: 


B  •  FS 


} 

ac 


(32) 


In  Eq.  (33),  N  is  the  number  of  teeth  on  the 
gear,  and  R^  is  the  raiiability  of  a  single  tooth 
on  the  gear.  Tha  reliability  of  any  tooth  on  a 
gear  is  equal  to  tha  raiiability  of  any  other  tooth 
on  the  gear.  To  transform  Eq.  (33)  into  a  life 
relationship,  substitute  Bq.  (28)  for  the  two  reli¬ 
abilities  into  the  log  of  Eq.  (33). 


‘  10, a 


i/b 


(34) 


Tha  gear  life,  I  ,  has  units  of  million 
gear  rotations. 


The  gear  dynamic  capacity,  C  ,  is  found  by 
substituting  Eq.  (34)  into  Eq.  (2^)  for  tha  tooth 
to  obtain  tha  analog  of  Bq.  (29)  for  the  gear. 
This  produces: 


C 
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Ct 


(35) 


Gear  Design 

Consider  the  design  of  a  gear  set  to  transmit 
80  hp  from  a  shaft  turning  at  5000  rpm  to  an  output 
shaft  turning  at  2500  rpm.  The  life  of  the  gears 
as  a  pair  in  hours  should  be  maximized  at  a  center 
distance  of  5.0  in. 

Table  1  is  the  program  output  listing  of  the 
problem  defining  constants  and  the  design  variable 
range  and  initial  value.  The  gears  are  to  have  a 
20*  nominal  pressure  angle  and  be  made  of  high 
strength,  heat  treated  steel  with  a  tooth  surface 
finish  of  32  rms.  The  material  surface  constant  of 
9800  psi  corresponds  to  a  surface  compression 
endurance  strength  of  200  000  psi  at  10^  fatigue 
cycles  and  a  raiiability  of  90  percent.  The  load- 
life  factor  of  8.93  is  from  the  ANSI/AGHA  2001  B88 
Standard,^*  and  tha  Weibull  slope  of  2.5  is  from 
the  NASA  Lewis  gear  test  data.*'^'  The  three 
design  variables  to  be  found  are  the  nutid>er  of 
teeth  on  the  pinion,  the  diametral  pitch  and  the 
face  width  of  the  gears.  For  the  gear,  the  number 
of  teeth  will  be  twice  that  of  tha  pinion. 

The  gears  should  have  a  face  width  which  gives 
the  pinion  a  length  to  pitch  diameter  ratio  between 
0.2  and  0.5.  Additional  strength  limits  placed  on 
the  design  are:  a  tooth  bending  stress  limit  of 
40  000  psi,  a  surface  contact  stress  limit  of 
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150  000  psl,  a  pressure  times  velocity  factor  of 
100  million  psi-ft/min,  and  an  oil  flash  temper¬ 
ature  limit  of  275  “f.  Involute  interference 
should  be  avoided  as  well.  These  limits  are  shown 
in  Table  2  which  lists  the  design  constraints  and 
merit  function  weighting  factors  in  the  output  data 
file. 

The  controlling  constraints  are  the  request 
for  a  center  distance  no  larger  than  5.0  in.  and 
the  maximum  pinion  length  to  diameter  ratio. 
Included  in  the  list  of  constraints  are  properties 
which  have  limits  that  any  design  will  satisfy, 
such  as:  pinion  torque,  transmitted  load,  dynamic 
load,  pinion  life,  and  the  oil  flash  temperature. 
These  constraints  enable  the  program  to  analyze 
each  design  for  their  values  and  display  them  as 
part  of  the  design  summary  for  each  optimum  and 
check  design.  A  velocity  factor  equation  calcu¬ 
lates  the  dynamic  load  as  prescribed  by  ANSI/AGMA 
2001  B88  with  a  quality  number  of  9.  All  calcu¬ 
lated  stresses  and  contact  pressures  use  the 
dynamic  load. 

Table  3  summarizes  the  program's  optimal  solu¬ 
tion  for  the  maximum  life  design  including  the 
design  variable  values,  the  merit  function  param¬ 
eter  values  and  all  constraint  values  and  limits. 
For  this  ideal  design,  the  mesh  life  is  28  000  hr. 
Table  4  gives  a  similar  summary  for  a  selected 
maximum  life  design.  Figure  2  shows  this  design 
which  has  46  and  92  teeth  on  the  two  gears,  a 
dieunetral  pitch  of  14  and  a  face  width  of  1.625  in. 
It  has  a  life  of  18  000  hr  with  all  other  con¬ 
straints  satisfied.  The  design  has  a  pinion  weight 
of  3.9  lb,  a  transmitted  load  of  650  lb  and  a  total 
dynamic  load  of  1840  lb. 

To  illustrate  the  flexibility  of  the  design 
program,  the  design  is  changed  to  request  a  minimum 
life  of  2000  hr  to  see  how  compact  the  design  can 
become. 

Six  changes  are  made  to  the  input  file.  The 
center  distance  limit  changes  from  an  upper  bound 
of  5.0  in.  to  a  lower  bound  of  0.0  in.  The  mesh 
life  limit  changes  from  a  lower  bound  of  0.0  hr  to 
a  lower  bound  of  2000  hr.  The  merit  function 
direction  changes  from  MAX  to  MIN.  And,  the  non¬ 
zero  weighting  factor  switches  from  mesh  life  to 
center  distance.  Labeling  changes  for  the  file 
itself  and  the  problem  title  make  up  the  last  two 
changes  which  enable  the  program  to  preserve  and 
distinguish  the  input  and  output  records  for  the 
two  design  requests. 

Table  5  summarizes  the  ideal  optimum  design 
for  the  minimum  size  design.  In  this  design,  the 
gear  center  distance  is  reduced  to  below  5.5  in.  by 
the  request  to  reduce  the  mesh  life  by  an  order  of 
magnitude.  Table  6  summarizes  the  nearest  practi¬ 
cal  design,  which  is  shown  in  Fig.  3.  In  this 
design,  the  pinion  and  the  gear  have  36  and  72 
teeth,  and  both  have  a  diametral  pitch  of  12  and  a 
face  width  of  1.5  in.  The  center  distance  is 
4.5  in.  and  the  mesh  life  is  2200  hr. 

Although  the  bending  and  contact  stresses  have 
increased  over  those  for  the  maximum  life  design, 
they  are  still  within  the  acceptable  criteria.  The 
scoring  limit  values  also  have  increased  to  accept¬ 
ably  larger  values.  Had  any  of  these  constraints 
been  reached,  the  bounding  constraints  would  have 
influenced  the  ideal  and  final  optimal  designs. 


It  should  be  noted  that  all  the  designs  have 
the  maximum  possible  face  width.  Designs  based  on 
pinion  weight  rather  than  center  distance  also 
required  the  maximum  possible  face  width.  Minimiz¬ 
ing  pinion  weight  for  a  life  of  2000  hr  produced 
the  same  optimal  design  as  that  found  by  minimizing 
the  gear  center  distance. 

Conclusions 

By  combining  the  power  of  optimization  with 
the  access  of  the  desk  top  computer,  a  practical 
spur  gear  life  design  program  has  been  written. 

The  program  uses  a  fixed  step,  modified  feasible 
directions  algorithm  to  search  for  the  optimum 
design.  Basic  relations  for  the  optimization 
algorithm  are  presented. 

Extensive  labels  and  keyboard  interactions 
give  the  designer  a  record  of  the  ideal  optimal 
design  and  any  user  specified  designs  entered  after 
the  optimum  has  been  found  and  reported.  Small 
changes  in  the  input  data  file  enable  one  to  redi¬ 
rect  the  design  objective  from  maximum  life  to 
minimum  size  and  to  change  the  controlling  design 
constraints . 

Two  gear  analysis  subroutines  and  the  compat¬ 
ible  input  data  file  apply  the  optimization  proce¬ 
dure  to  the  gear  life  design  problem.  The  first 
subroutine  analyzes  a  gear  design  for  the  con¬ 
straint  variable  values.  And  the  second  subroutine 
analyzes  a  gear  design  for  the  merit  function 
parameter  values.  The  input  data  file  defines 
these  values,  sets  the  limits  for  the  design  con¬ 
straints,  and  selects  the  active  merit  function 
parameters. 

Tooth  surface  pitting  fatigue  life  produces 
the  finite  life  of  the  gear  set.  A  fatigue  life 
model  for  this  mode  of  failure  is  presented.  Sta¬ 
tistical  variations  in  gear  life  as  predicted  by 
the  two  parameter  Weibull  distribution  are  given 
also. 

Two  designs  illustrate  the  design  calculations 
of  the  program.  One  gear  set  is  designed  for  maxi¬ 
mum  life  with  a  required  center  distance,  and  a 
second  design  is  obtained  for  a  minimum  center 
distance  at  a  minimum  acceptable  life.  The  design 
procedure  considers  other  constraints  as  well  to 
obtain  practical  gear  designs.  Minimum  weight  and 
minimum  center  distance  designs  were  found  to  be 
identical . 
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TABLE  1.  -  GEAR  DESIGN  CONSTANTS  AND  VARIABLES 
FOR  SPUR  GEAR  LIFE  MAXIMIZATION 

[Design  with  modified  gradient  optimization  using  a  maximum  step 
limit  and  scaled  variables.] 

_  (a)  Fixed  design  requirements 


Poisson's  ratio  .  0.25 

Elastic  modulus,  psi .  3x10^ 

Pressure  angle,  deg  .  20 

Gear  ratio .  2 

Transmitted  power,  hp  .  80 

Pinion  speed,  rpm  .  5000 

Material  weight  density,  Ib/in.^  .  0.283 

Material  surface  constant,  psi  .  9800 

Weibull  slope  .  2.5 

Load-life  factor  .  8.93 

Reliability  .  0.9 

Base  temperature,  ®F .  120 

Tooth  surface  fiish,  rms .  32 


(b)  Estimated  values  of  the  three  independent 
design  variables 


Low 

High 

Initial 

Pinion  teeth 

10 

100 

40 

Di^unetral  pitch,  in.*^ 

4 

28 

14 

Face  width,  in. 

0.5 

5 

2.5 
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TABLE  2 


16  MAXIMUM  LIFE  CONSTRAINTS  AND  MERIT  FUNCTION 


(a>  The  16  constraint  limits 


Constraint 

Limit 

Type  of 

bound 

Involute  interference,  in. 

0.001 

Lower 

Lower  face  width-to-diameter  ratio 

0.2 

Lower 

Upper  face  width-to-diameter  ratio 

0.5 

Upper 

Pinion  weight,  lb 

0 

Lower 

Center  distance,  in. 

5.0 

Upper 

Pinion  torque,  lb-in. 

0 

Lower 

Transmitted  load,  lb 

0 

Lower 

Total  dynamic  load,  lb 

0 

Ixjwer 

AGMA  bending  stress,  psi 

0.4x10® 

Upper 

Full  load  contact  stress,  psi 

1.5x10® 

Upper 

Gear  tip  hertz  pressure,  psi 

1.5x10® 

Upper 

Pinion  life,  cycles 

0 

Lower 

Mesh  life,  hr 

0 

Lower 

Pitch  line  velocity,  ft/min 

0 

Lower 

PV  factor,  psi-ft/min 

100 

Upper 

Flash  temperature,  *F 

275 

Upper 

(b)  Maximize  the  objective  function  (OBJ),  which 


has  three  terms.  OBJ  =■  the  linear  sum  of 


Component 

Unit 

Multiplied 

by 

Mesh  life 

in  thousands  of  hours 

1 

Center  distance 

in  inches 

0 

Pinion  weight 

in  pounds 

0 

TABLE  3.  -  MAXIMUM  LIFE  OPTIMAL  DESIGN 
(Optimization  successful  in  50  steps.) 

(a)  The  final  design  vectors 


Design  parameter 


Pinion  teeth  -0.24857  43.81443 
Diametral  pitch,  in."'  -.23797  13.14433 
Face  width,  in.  -.48148  1.66667 


(b)  Components  of  the  maximum  objective 
function,  28.1361 


Value  Multiplied 
by 


(c)  The  16  constraint  values 


Component 

Value 

Mesh  life,  hrxlO® 
Center  distance,  in. 
Pinion  weight,  lb  1 

28.136 

5.0000 

4.1161 

Constant 

Value 

Limit 

Type  of 
bound 

Involute  interference,  in. 

0.23224 

1x10"® 

Lower 

Lower  face  width-to-diameter  ratio 

0.5 

0.2 

Lower 

Upper  face  width-to-diameter  ratio 

0.5 

0.5 

Upper 

Pinion  weight,  lb 

4.1161 

0 

Lower 

Center  distance,  in. 

5.0 

5.0 

Upper 

Pinion  torque,  lb-in. 

1008.4 

0 

Lower 

Transmitted  load,  lb 

643.87 

0 

Lower 

Total  dynamic  load,  lb 

1825.9 

0 

Lower 

AGMA  bending  stress,  psi 

0.34422x10® 

0.4x10® 

Upper 

Full  load  contact  stress,  psi 

1.2290x10® 

1.5x10® 

Upper 

Gear  tip  hertz  pressure,  psi 

0.80590x10® 

1.5x10® 

Upper 

Pinion  life,  cycles 

9.5274x10’ 

0 

Lower 

Mesh  life,  hr 

28.136x10® 

0 

Lower 

Pitch  line  velocity,  ft/min 

4100.2 

0 

Lower 

PV  factor,  psi-ft/min 

67.268x10® 

100x10® 

Upper 

Flash  temperature,  *F 

203.96 

275 

Upper 
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TABLE  4.  -  MAXIMUM  LIFE  DESIGN  CHECK 


(a)  Daalgn  check 


Design  parameter 


Pinion  teeth 
Diametral  pitch,  in. 
Face  width,  in. 


(b)  The  components  of  the  maximum  objective 
function,  18.1249 


46.000 

14.000 

1.625 


Component 


Mash  life,  hrxio’ 
Center  distance,  in. 
Pinion  weight,  lb 
(c)  The  16  constraint  values 


TABLE  S.  -  MINIMUM  SIZE  OPTIMAL  DESIGN 
(Optimization  successful  in  27  steps.) 

(a)  The  final  design  vectors 


Design  parameter 


Pinion  teeth  -0.45045 
Diametral  pitch,  in.”^  -.36590 
Face  width,  in,  -.55778 


(b)  The  components  of  the  minimum  objective 
function,  4.48739 


34.72992 

11.60917 

1.49501 


Component 

Value 

Mesh  life,  hrxio’ 
Center  distance,  in. 
Pinion  weight,  lb 

2.0236 

4.4874 

2.9739 

Multiplied 

by 


(c)  The  16  constraint  values 


Value  Multiplied 
by 


18.125 

4.9286 

3.8993 


Constraint 

Value 

Limit 

Type  of 
bound 

Involute  interference,  in. 

0.23810 

1x10'* 

Lower 

Lower  face  width-to-diamater  ratio 

0.49457 

0.2 

Lower 

Upper  face  width-to-diameter  ratio 

0.49457 

0.5 

Upper 

Pinion  weight,  lb 

3.8993 

0 

Lower 

Center  distance,  in. 

4.9286 

5.0 

Upper 

Pinion  torque,  lb-in. 

1008.4 

0 

Lower 

Transmitted  load,  lb 

653.20 

0 

Lower 

Total  dynamic  load,  lb 

1843.8 

0 

Lower 

AGMA  bending  stress,  psi 

0.37648xl0’ 

0.4x10* 

Upper 

Full  load  contact  stress,  psi 

1.2582x10* 

1.5x10* 

Upper 

Gear  tip  hertz  pressure,  psi 

0.81864x10* 

1.5x10* 

Upper 

Pinion  life,  cycles 

6.1374x10® 

0 

Lower 

Mesh  life,  hr 

18.125x10* 

0 

Lower 

Pitch  line  velocity,  ft/min 

4041.6 

0 

Lower 

PV,  factor,  psi-ft/min 

64.271x10* 

100x10* 

Upper 

Flash  temperature,  ”F 

202.37 

275 

Upper 

Constraint 

Value 

Limit 

Type  of 
bound 

Involute  interference,  in. 

0.20055 

1x10’* 

Lower 

Lower  face  width-to-diameter  ratio 

0.49919 

0.2 

Lower 

Upper  face  width-to-diameter  ratio 

0.49919 

0.5 

Upper 

Pinion  weight,  lb 

2.9699 

0 

Lower 

Center  distance,  in. 

4.4870 

0 

Lower 

Pinion  torque,  lb-in. 

1008.4 

0 

Lower 

Transmitted  load,  lb 

717.48 

0 

Lower 

Total  dynamic  load,  lb 

1965.3 

0 

Lower 

AGMA  bending  stress,  psi 

0.38256x10* 

0.4x10* 

Upper 

Full  load  contact  stress,  psi 

1.4325x10* 

1.5x10* 

Upper 

Gear  tip  hertz  pressure,  psi 

0.98539x10* 

1.5x10* 

Upper 

Pinion  life,  cycles 

0.67728x10® 

0 

Lower 

Mesh  life,  hr 

2.0001x10* 

2.0x10* 

Lower 

Pitch  line  velocity,  ft/min 

3679.5 

0 

Lower 

PV  factor,  psi-ft/min 

92.242x10* 

100x10* 

Upper 

Flash  temperature,  *F 

234.28 

275 

Upper 
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TABLE  6.  -  MINIMUM  SIZE  DESIGN  CHECK 


(a)  The  final  design  vectors 


Design  parameter 


Pinion  teeth 

36.000 

Diametral  pitch,  in.“* 

12.000 

Face  width,  in. 

1.500 

(b)  The  components  of  minimum  objective 
function,  4.500 


Constraint 

Value 

Limit 

Type  of 
bound 

Involute  interference/  in. 

0.20588 

1x10'* 

Lower 

Lower  face  width-to*diameter  ratio 

0.5 

0.2 

Lower 

Upper  face  width-to-diameter  ratio 

0.5 

0.5 

Upper 

Pinion  weight,  lb 

3.0006 

0 

Lower 

Center  distance,  in. 

4.5 

0 

Lower 

Pinion  torque,  lb-in. 

1008.4 

0 

Lower 

Transmitted  load,  lb 

715.41 

0 

Lower 

Total  dynamic  load/  lb 

1961.4 

0 

Lower 

AGMA  bending  stress,  psi 

0.38980x10^ 

0.4x10* 

Upper 

Full  load  contact  stress,  psi 

1.4239x10* 

1.5x10* 

Upper 

Gear  tip  hertz  pressure,  psi 

0.97093x10* 

1.5x10* 

Upper 

Pinion  life,  cycles 

0.74373x10’ 

0 

Lower 

Mesh  life,  hr 

2.1964x10* 

2.0x10* 

Lower 

Pitch  line  velocity,  ft/min 

3690.2 

0 

Lower 

PV  factor,  psi-ft/min 

88.071x10* 

100x10* 

Upper 

Flash  temperature,  "f 

230.41 

275 

Upper 
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